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We study the exact solution for two atomic particles in an optical lattice interacting via a Fes- 
hbach resonance. The analysis includes the influence of all higher bands, as well as the proper 
renormalization of molecular energy in the closed channel. Using an expansion in Bloch waves, we 
show that the problem reduces to a simple matrix equation, which can be solved numerically very 
efficient. This exact solution allows for the precise determination of the parameters in the Hub- 
bard model and the two-particle bound state energy. We identify the regime, where a single band 
Hubbard model fails to describe the scattering of the atoms as well as the bound states. 



Cold atomic gases in optical lattices represent a perfect 
laboratory system for the quantum simulation of strongly 
correlated many-body systems described by Hubbard 
models [H, 0, 0] . Recently, experimental and theoretical 
efforts focus on the observation of a Fermionic Mott insu- 
lator 0, [B| , and the ultimate goal towards the realization 
of magnetic and superconducting phases. The quanti- 
tative understanding of these experimental results and 
the comparison with the theoretical predictions require 
a precise knowledge of the parameters in the Hubbard 
model for cold atomic gases interacting with a Feshbach 
resonance @. In this letter, we present the solution to 
the two-particle problem in an optical lattice interact- 
ing via a Feshbach resonance, and provide a microscopic 
derivation of the parameters in the Hubbard model and 
the two-particle bound state energies. 

The two-particle interaction potential between parti- 
cles in ultra-cold atomic gases is well described by the 
pseudo-potential or in the presence of a Feshbach res- 
onances within a two-channel model 0, 0]. The two- 
particle problem within confined geometries has exten- 
sively been studied for the one-dimensional setup with 
strong transverse confining 0, 0] , and the harmonic trap- 
ping potential [Io| . In addition, the influence of optical 
lattices has been studied for the deep lattices, where the 
influence of higher bands have been included semiclassi- 
caly ll| or using the exact solution for a harmonic oscil- 



lator within each well of the lattice [12|, |_13|, Here, 
we analyze the two-particle problem interacting via a 
Feshbach resonance in a three-dimensional optical lat- 
tice and show that the equations can be efficiently solved 
numerically. The solution provides the exact scattering 
properties and bound state energies of two-particles in 
an optical lattice of arbitrary strength, see Fig. [TJ 

From the exact scattering amplitude, we find the mi- 
croscopic derivation for the interaction parameters in the 
Hubbard model. The simplest Hubbard model describes 
bosonic particles with creation (annihilation) operators 
b\ (bi), and on-site interaction U (extension to fermionic 
particles with spin is straightforward), 
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FIG. 1: Exact bound state energies (red line) for two-atoms 
in an optical lattice of strength V = 12E r and K = 0. The 
additional bound states (dashed grey lines) are weakly cou- 



pled to atoms in the lowest Bloch band, i.e. 



< to- 
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The green dotted line denotes the bound state energies pre- 
dicted from the Hubbard model with the on-site interaction 
determined by Eq. (TTU)) ; its deviations from the exact bound 
state indicate the break down of the Hubbard model. The 
blue dotted line denotes the bound states neglecting the cor- 
rection xo — |^o| 2 G(0) = —3.5. 



The hopping energies iy derive from a single parti- 
cle band structure calculation, and are related to the 
dispersion relation in the lowest Bloch band E^q) — 
— cosq(Ri — Rj). In turn, the on-site interac- 
tion U is conventionally derived for weak interaction 
strengths and deep optical lattices by replacing the exact 
pseudo-potential by a ^-function interaction and restrict- 
ing the system to the lowest Bloch band [l5(; the latter 
step corresponds to introducing a short distance cut-off 
A comparable to the lattice spacing a. This approach is 
restricted to weak interactions with |a s |/a <C 1; here, a s 
is the s-wave scattering length. In the general situation, 
the precise derivation of the interaction potential in the 
Hubbard model for arbitrary interaction is obtained by 
comparing the exact scattering properties of two parti- 
cles in an optical lattice with the scattering amplitude 
predicted from the Hubbard model. This approach is 
in analogy to the description of the interaction bewteen 
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cold gases in free space in terms of a pseudo-potential: 
the strength of the pseudo-potential is fixed by the con- 
dition to reproduce the exact scattering properties. 

In the following, the interaction between the two- 
particles is given by a Feshbach resonance, which can 
be conveniently described by the two-channel approach. 
Then, the Feshbach resonance is characterized by the de- 
tuning v and the coupling g between the open and closed 
channel, and gives rise to the scattering amplitude [3] 
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with k the incoming momentum and fi the reduced mass. 
The scattering length a s = —2ng 2 /A'Kh 2 v and the effec- 
tive range r b = nh A / n 2 g 2 are experimentally accessible 
by measuring the bound state energy of the molecules 
across the Feshbach resonance [l6| ■ 

The two particles in the open channel are described by 
the wave function ^(x, y) with x and y the position of 
the particles. In order to capture the above characteris- 
tics of a Feshbach resonance, it is enough to describe the 
closed channel by a single molecular state </>(z). Then, 
the Schrodinger equation for the energy eigenstates re- 
duces to 

[E - HS + H b Q ] ^(x, y) = g J <iza(r)0(z)<5 (z - R) (3) 
[E-v Q -H™]4>{z) = g J d X dya(r)V(x,y)5( Z -R), 

where the single particle physics is described by the 
Hamiltonians Hq = —^-A + V^(x) with V a the optical 
lattices and the molecular mass m m = m a + mi,. Fur- 
thermore, we have introduced the relative r = x — y and 
center of mass coordinates R = (ra a x + m;,y) / (m„ + . 
The properties of the Feshbach resonance are determined 
by the coupling strength g and the bare detuning vo, 
while a(r) = exp(-r 2 /2A 2 )/(2A 2 7r) 3 / 2 — > S(r) accounts 
for a regularization of the coupling with cut-off A. In 
the limit A — > 0, the bare detuning Uq entering the 
microscopic theory is related to the physical observable 
detuning v via v a = v — u Ten with the renormalization 
v Icn = - 9 V(2fr 2 7r 3 / 2 A). 

The periodic structure of the optical lattice is charac- 
terized by the lattice vectors {Rj}- The single particle 
properties are then fully determined by the Bloch wave 
functions V'^k a ( x )' ^m.qiy)^ an d 0s.k(z) with the cor- 
responding band energies E^(k a ), E^(kb), and £™(K); 
here, k , kf, and K are the quasi-momentum, while n, 
m, and s characterize the different Bloch bands. In 
the following, we measure energies with respect to the 
ground state energy of two-particles in the lowest Bloch 
band, i.e., E§(0) + £$(0) = 0. The discrete translation 
invariance provides the conservation of the total quasi- 
momentum K = k a + kf,. Then, the general solution 



with fixed total quasi-momentum K can be written as 
^(x,y) = ^=^^^" l (q)V^ q (x)^ K _ q (y), 



N 

and 4>(z) = X] s ^ s< ^s,k(z). Inserting this expansion in 
Eq. ([3]), we obtain 

[£-£„ m (q)]^ m (q) = ^]>>r(q)^ (4) 

s 

[E-v - E? (K)] R s = w ]T 1 ]T K m (aW im (q) ■ 

n,m q 

Here, we have introduced the notation E nm (q) = 
^n(q) + ^m(K — q), and the characteristic coupling en- 
ergy w — g/^/Vo with Vq the volume of the unit cell. The 
dimensionless coupling elements reduce to 



/i" m (q) 
VNVc 



U = J dxdy [C, q (x)< K - q (y)] *«(r)0 s , K (R) , 



with the notation h s nm (q) = [h" m (q)}*, while NV de- 
notes the quantization volume. Substituting the bare 
detuning vq with the physical detuning v by adding on 
both sides of Eq. ([4J the renormalization z/ ron , we obtain 



\E-v- ET] R s ~w 2 J2 X S t(E)R f ' = 0. 



(5) 



The matrix \t describes the shift of the Feshbach res- 
onance due to the change in dispersion relation of the 
particles in the open channel; this phenomena is in anal- 
ogy to the lamb shift of atoms in a cavity It takes 
the form (vq denotes the volume of the Brioullinc zone) 



^ m (q)/ir(q) , K m (q)hr(<l) 



E - E nm {q) + i<q 



E nm {q) 



The quantities E nm (q) and ft.™ m (q) are the energies and 
coupling parameters for the system in absence of an opti- 
cal lattice. The first term in the above equation describes 
the influence of higher bands, while the second term ap- 
pears from the renormalization. The divergent parts in 
the two terms cancel each other, and \t remains finite 
in the limit A — > 0. This behavior can be easily under- 
stood: for large Bloch bands, the influence of the optical 
lattices vanishes and the coupling elements /i™ m (q) and 
energies E nm (q) reduce to the values of the free system. 
Then, the terms in the bracket cancel each other, and the 
summation over the Bloch bands converges. For a finite 
short distance cut-off A, the corrections vanishes with 
~ A; i.e., the convergence is very slow in the number of 
Bloch bands. 

In the following, we discuss the setup with a three- 
dimensional cubic lattice Vo-(x) = V a Y^=i sin 2 (fc L 2!i) 
with lattice spacing a — 7r/fc L and recoil energy E r = 
h 2 k 2 /2m. For equal particle species m a = nib an d far 
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detuned optical lattice, the relative strengths of the lat- 
tice potentials naturally satisfy V a = Vb = V m /2, We fo- 
cus on a wide Feshbach resonance; the generalization to 
a narrow Feshbach resonance is straightforward. A wide 
Feshbach resonance is obtained in the limit v, g — > oo 
with a fixed s-wave scattering length a s = —mg 2 /4irh 2 is, 
and the energies E and E™ in the first term in Eq. ([5]) 
can be dropped. 
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FIG. 2: Bound state energies (red line) for different lattice 
depths: (a) For weak optical lattices with V = 4_E r , the 
Bloch bands in three-dimensions nearly overlap, (b) For deep 
optical lattices with V = 50i5 r , the bound states energies 
are compared with the energies obtained by replacing the 
optical lattice by a harmonic well with trapping frequency 
uj p — ^/4VE r /h (dashed blue line). 

Bound states: The equation for the energies E B of 
the repulsive and attractive bound states reduces to the 
eigenvalue equation 



E 
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■5t + W X 8 t{E B )]R t = 
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with R s the molecular wave function of the bound state, 
and W — —g 2 /vVa = (8/n)E r a s /a; the numerical so- 
lution is shown in Fig. Q] and Fig. [2] For a fixed value 
of the s-wave scattering length a s and quasi-momentum 
K, there appear now several bound states. This behav- 
ior is in strong contrast to the free system, where for a 
fixed center of mass motion only a single bound state 
exists for a repulsive scattering length a s . Here, the ap- 
pearance of several bound states is a consequence of the 
reduced translation symmetry: molecular states and two- 
particles states in the open channel differing in center of 
mass motion by a reciprocal lattice vector are coupled 
the periodic lattice. These coupling strength are given 
by the overlaps /i" m (q). 

The matrix elements \t (E) are fully determined by the 
single particle properties such as the Bloch wave func- 
tions and Band structure, and can be efficiently deter- 
mined numerically. For the three-dimensional cubic lat- 
tice above, the single particle wave function separate for 
each space direction, i.e., u n> k(x) = Yl i=1 u^. f. ( x i)> an d 
it is therefore sufficient to determine the Bloch wave func- 
tions and energies for a one-dimensinal setup; the Bloch 
wave functions are determined using Z — 20 1 3 recipro- 
cal lattice vectors. The integration and the summation 
over the different Bloch bands is performed. While the 



integration converges very quickly using N = 41 3 unit 
cells, the limiting factor in accuracy for the matrix ele- 
ments \t O^O ls the slow convergence with the number of 
Bloch bands: the restriction to a finite number of Bloch 
bands in the summation corresponds to introducing a 
high energy cut-off A ~ l/(fc L S' 1 / 3 ). Consequently, the 
summation converges with ~ l/.? 1 / 3 , and a finite size 
scaling analysis can be performed; the number of Bloch 
bands included in this analysis was S — ll 3 . Then, the 
matrix elements Xt can be calculated with an accuracy 
better than 1%; the convergence has been extensively 
tested for varying number of unit cells N = 41 3 . . . 201 3 , 
S = ll 3 . . . 2I 3 , and Z = 100 3 . . . 500 3 . 

Scattering amplitude: Next, we analyze the scattering 
states in the lowest Bloch band; 



^ nm (q) = vr(q) + 



A" m (q, k,K) 



E-E nm (q,K) 



(T) 



with <po(q) = <^n,0^ro,o<5q,k an incoming wave at rela- 
tive momentum q and center of mass momentum K in 
the lowest Bloch band. Then the generalization of the 
s-wave scattering length in free space is obtained via 
A = A 00 (q — > 0, 0, 0) at low energy E — > of the incoming 
wave. Again, the scattering amplitude A is fully deter- 
mined by the matrix xl ■ Introducing the notation R s a for 
the eigenvectors of the matrix xi with eigenvalues Xa , the 
scattering amplitude reduces to (W = (8/ir)E r a s /a) 



W\w \ 



W Xa 1 - W X0 



(8) 



with w a = ^2 S /i°°(0)i?* the width of each scattering res- 
onance. The crossing of each bound state with the lowest 
Bloch band, see Fig. [H gives rise to a pole in the scat- 
tering amplitude and describes a scattering resonance. 
Except for the first resonance, the couplings are in gen- 
eral weak and the scattering amplitude is dominated by 
the lowest eigenvalue xo & n d width wo. However, these 
additional resonances can give rise to characteristic loss 
features for cold atoms in an optical lattice at large s- 
wave interactions. 

In the following, we will now compare this exact value 
for the scattering amplitude for two particles in an opti- 
cal lattice with the predictions from the Hubbard model 
Eq. JTJ). For an on-site interaction U the scattering solu- 
tion in the Hubbard model takes the form 1181 



^HM(q) 



A, 



3q,k 



£-£ 00 (q,K) 



irj 



(9) 



where A HM describes the scattering amplitude in the 
Hubbard model; A HM = U/[l - UG{E)] with G(E) = 
I [E ~ Eoo(q) + ii]}~ 1 ■ F° r nearest neighbor hopping 
t and low scattering energies E — ► 0, G(0) reduces to 
G(0) = c/2t with c w -0.2527. The effective on-site 
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0.0855 


-4.188 


2.412 


-0.6/ E r 
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0.0308 


-26.82 


5.954 


-2.3/ E r 


12 


0.0122 


-101.3 


9.483 


-3.5/ E r 


16 


0.00533 


-303.0 


12.63 


-3.7/E r 


20 


0.00249 


-788.2 


15.50 


-1.8/ E r 



TABLE I: Effective parameter in the Hubbard model with 
hopping t. The on-site interaction U is given by Eq. (1101) . 



interaction U is therefore completely fixed by the condi- 
tion, that the Hubbard model reproduces the exact two- 
particle properties, i.e., A HM = A, and we obtain 



U 



1 



W\w \ 



x- 



G(0) l-W[ X o-\w \ 2 G(p)Y 



(10) 



The parameters for different strengths of the optical lat- 
tice are shown in Table Q] The contribution VF|u>o| 2 de- 
scribes the dominant part for weak interactions, while the 
correction \o — \wq\ 2 G(0) becomes relevant for stronger 
interactions. It is important to stress, that this deriva- 
tion of the on-site interaction U is valid for arbitrary 
values of the s-wave scattering length a s , and gives rise 
to a finite value for a s — > ±oo. However, its va- 
lidity is restricted to low scattering energies: first, ad- 
ditional interaction terms beyond the on-site interaction 
U can play an important role and will account for the 
full momentum dependence of the scattering amplitude 
A" m (q,k,K). Second, the bound state energies can be 
strongly modified by additional terms, which are not in- 
cluded in a single band Hubbard model. A test for the 
validity of the Hubbard model is therefore the compar- 
ison with the exact bound state energy and the repul- 
sive/attractive bound states predicted from the Hubbard 
model. The bound states within the Hubbard model are 
determined by poles in the scattering amplitude A HM , ke., 
UG(E) — I. A comparison with the exact bound state 
energies is shown in Fig. [U and we find already very 
strong deviations for a s /a > 0.02 at V = 12E r : the va- 
lidity for the description of bound states in the Hubbard 
model is limited to very weak interactions. 

Finally, for deep optical lattices V/E r > 1, the width of 
the lowest Bloch band provides a small parameter char- 
acterized by the the hopping energy i/B r C 1. As a 
consequence, for all energies E of the order of the band 
width E r*> 12t, the first term in the matrix xt dominates 



rfq 



v E- E oa (q) +ir) 
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while all the remaining terms from higher Bloch bands 
as well as the renormalization provide a contribution 
~ 1/E r . Consequently, the results reduces to the well 
known approach [h| for the derivation of the Hubbard 
parameters, where the influence of higher bands are ne- 
glected and the pseudo potential is replaced by a 5- 



function. Then, the momentum dependence of the in- 
teraction potential in the Hubbard model reduces to 



£/(q,k,K) = 
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{'Po, K- 



>0. 



This term also accounts for contributions such as nearest- 
neighbor interactions and correlated hopping [HI, [l|| [2(| • 
For increasing interactions, the shift xo~ |wo| 2 G(0) ~ 
1/E r in Eq. pop becomes important, and in addition the 
bound state energies E B start to deviate from the predic- 
tions within the Hubbard model. The crossover from the 
two regimes can be self-consistently checked: the higher 
bands become relevant as soon as E B becomes in the 
range of the Bloch band separation U ~ hw p , while the 
renormalization requires W <C E ri i.e., a s <C a. In the 
limit of deep optical lattices, the first condition is always 
more stringent and reduces to a s /a <§; (Er/V) 1 ^ 4 /2y/n, 
which has previously been suggested [l5| . In order to de- 
rive Hubbard models which reproduce the bound state as 
well as the scattering states, the influence of the higher 
bands as well as the renormalization have to be included; 
in contrast to recent attempts to include the influence of 
higher bands alone [20]. 
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